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Abstract 

The superfluidity of neutron matter in the channel 1 5'o is studied by taking into account the 
effect of the ground-state correlations in the self-energy To this purpose the gap equation has 
been solved within the generalized Gorkov approach. A sizeable suppression of the energy gap 
is driven by the quasi-particle strength around the Fermi surface. 
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1. Introduction. 



It is widely recognized that superfluidity is an extremely subtle process when it is considered 
on an entirely microscopic level. Most of the calculations in nuclear physics, at least concerning 
finite nuclei, are therefore based on phenomenological effective interactions. On the other hand 
a series of nuclear matter and neutron matter calculations exists with the bare NN force used 
in the gap equation . This is insofar the first step of a systematic microscopic approach as it 
is well known that in the gap equation the bare particle-particle interaction has to be taken to 
the lowest order and not a G-matrix, i.e. a ladder summation, since the gap equation is in itself 
a two-body equation [0, pi. The next terms in the effective pairing interaction are screening 
terms. These are due to the possibility of medium polarization which, when treated within the 
induced interaction theory [|[], has turned out to strongly affect the pairing gap ]5], [/J. On 
the other hand the self-energy corrections due to polarization effects have attracted much less 
attention. Only very recently it has been stressed that the consideration of the quasi-particle 
strength can have a sizeable effect on the gap value ||. Here we want to extend the latter work 
in investigating correlation effects beyond the ones taken into account in Ref. || and we will see 
that they are indeed quite important. Our study must be considered as an intermediate step 
towards a fully consistent treatment where self-energy and vertex corrections are taken into 
account on an absolutely equal footing. The latter aspect may turn out to be decisive because 
sometimes strong cancellation between the two contributions can occur (see e.g. discussion in 
Ref. 0). All those considerations are of great importance for a more microscopic understanding 
of superfluidity in nuclei as well as in neutron stars, but also for a precise estimate of the pairing 
gap. In fact, pairing in exotic nuclei at present is studied with effective density dependent 
interactions modeled on pairing calculations in neutron matter P, ETUI. In the latter systems it 



is well known that superfluidity drives their rotational dynamics [|TTJ as well as their cooling 



12]. In this work we shall study pure neutron matter but we expect that very analogous effects 



will occur in symmetric nuclear matter. 
2. Neutron self-energy. 

Let us first discuss some properties of the single-particle self-energy T, p (u) of neutron matter. 



In the Brueckner approach |TB|] the perturbative expansion of E can be recast according to the 
number of hole lines as follows 

E p (w) = £» + E» + . . . (1) 

The on-shell values of S 1 represent the Brueckner-Hartree-Fock (BHF) mean field and the 
on-shell values of S 2 give the so called rearrangement term, which gives largest contribution to 
the ground-state correlations. The off-shell values enter several physical properties of neutron 
matter, including pairing. In terms of the self-energy one may calculate, at a given order of the 
hole-line expansion, the quasi-particle energy, as solution uj v of the equation 

p 2 

where is the Fermi energy. We here neglect the imaginary part of the self-energy. The 
quasi-particle energy around the Fermi surface is obtained expanding the self-energy around 
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p = Pf and to = 0: 



^ = p 2 -p\ = p 2 -p% m 
2m* 2m m P m n 



where m*,m e and m p are the effective mass, e-mass and p-mass, respectively. The latter two 
masses are defined as 

i - m 



m P = m 



m p = m 



.dcu 

. ( m /<9£\ ' 
Pf \dp J f. 



(4) 



(5) 



The partial derivatives are evaluated at the Fermi surface. The fc-mass is related to the non- 
locality of the mean field and, if the self-energy is ^-independent (static limit), it coincides 
with the effective mass. This quantity is of great interest whenever the momentum dependence 
of the mean field can give some effects such as transverse flows in heavy-ion collisions |fU 



The e-mass is related to the quasi-particle strength. This latter gives the discontinuity of the 
momentum distribution at the Fermi surface, and measures the amount of correlations included 
in the considered approximation. 

The definition of m*, fc-mass and e-mass can be extended to p ^ pf replacing in Eqs. @) 
and fl5|) pf by p and uj = by uj p . Their properties have been extensively studied in Ref. Hl3| . 

The self-energy has been calculated in two approximations: up to first order (BHF) and 
up to the second order (EBHF) of the hole-line expansion in the framework of the Brueckner 



theory adopting the continuous choice iT5| . The Argonne V14 potential has been used for the 
bare interaction [16]. The calculations have been performed for a range of Fermi momenta 
where the energy gap is expected to be the largest, i.e. 0.5 < k F < 1.3/m -1 , corresponding 
to a density range from 0.0042 through 0.074 fm~ 3 . A typical result for the off-shell neutron 
self-energy T, p {oj) is plotted in Fig. 1. The contribution from the rearrangement term shows 
a pronounced enhancement in the vicinity of the Fermi energy, which is to be traced back to 
the high probability amplitude for particle- hole excitations near At high momenta this 

contribution vanishes. From H p (uS) the effective masses are extracted according to Eq.s (4) and 
(5). They are depicted in Fig. 2, where the full calculation is compared to that including only 
the BHF self-energy. We may distinguish two momentum intervals: at k ~ kp the momentum 
dependence of the effective mass m* is characterized by a bump, whose peak value exceeds 
the value of the bare mass; far above kp the bare mass limit is approached. One should take 
into account that in this range of kp the neutron density is quite small (at the maximum 
kp = l.3fm~ 1 p = 0.074/m~ 3 ). This behaviour of the effective mass m* is due mostly to the 
e-mass, as shown in the lower panel of Fig. 2. In both panels of Fig. 2 it is also reported, for 
comparison, the effective mass in the BHF limit (only E 1 included), which exhibits a much 
less pronounced bump at the Fermi energy. It is precisely the influence of this increased bump 
structure, which we want to investigate here, since in Ref. || the rearrangement term has not 
been taken into account. 

3. Generalized gap equation 

The generalized BCS theory can be found in various textbooks on the many-body problem 



17, 18]. Here we follow closely the formalism developed in [H, where the gap equation is 
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written as 



A p (oj) 



The kernel T is defined as 



(6) 



r p (uj) 



G p (-u;)G» = [G;\u,)G;\-u>) + ^u,)]- 1 



(7) 



The functions G p (lv) and G p (u) are the nucleon propagators of neutron matter in the normal 
state and in the superfluid state, respectively . The ^-symmetry in the two propagators is to 
be traced back to the time-reversal invariance of the Cooper pairs. The effective interaction 
V is the block of all irreducible diagrams of the interaction. The short-range correlations 
(ladder diagrams) are already taken into account by the gap equation and do not appear in 
the irreducible block V. The long-range components have been studied in the context of the 
induced interaction approach [7J. In the present note we only consider the first term in the 
perturbative expansion of V, namely the bare interaction, because we want to disentangle the 
influence of correlations coming only from the self-energy expansion from the ones due to the 
induced interaction. The complete solution of the generalized gap equation requires a further 
effort . 

Assuming the pairing interaction to be identified with the bare interaction V p , p >, the energy 
gap does not depend on the energy (static limit), i.e. A p (u) = A p . In this limit the self-energy 
corrections are not expected to modify the analytical structure of the kernel r p (u;) which now 
is an even function of energy: at each momentum p there exist two symmetric poles ±Q p in 
the complex a;-plane. The u integration can be performed as follows: 



dui 
2k~i 



z 2 

V 



where we denote by Z 2 the residue of the kernel at the pole fl p . Since the largest contribution 
to the integral is coming from the pole part of the two Green's function, we expand the single 
particle propagator to first order in oj p (see Eq. (||)), G p 
the denominator in Eq. (|7j) becomes 



Z rt 



[uo — UJ r , 



and therefore 



3TV 



w2) + A 



(9) 



where Z 2 is given by 



and 



z: 
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In this approximation the generalized gap equation, Eq. (6), becomes 

1 r d 3 p' ZpVpp'Zp/ 



A, 



(27T) 



A p , 



(10) 
(11) 

(12) 
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where A p = Z P A P is the real pairing correction to the quasi-particle energy spectrum. The 
main difference from the BCS limit is the presence of the quasi-particle strength, which is less 
than one in a small region around the Fermi surface as we saw in the previous section. The 
pairing interaction turns out to be reduced in that region, where the Cooper pairs are mainly 
formed. This is the way self-energy corrections come into play suppressing the pairing gap. 
As to the self-energy effects, Eq. ( |T2"D is quite general because there is no expansion in the 
self-energy so far. The self-energy corrections can be taken at any order of approximation. 

One may further restrict the ^-integration to only the pole part at the Fermi energy, i.e. ex- 
panding the self-energy near the Fermi surface according to Eq. (3). In this case the integration 
can also be performed analytically and one easily obtains 

r 2 f d3 p' Vp tP iA p , 



A = _z 2 / —J- p ' p p (13) 

~' i 2 *) 3 ~ p F ) 2 /m* 2 + A 2 p ,' 



where Zp is the quasi-particle strength at the Fermi surface and coincides with the inverse of 
the e-mass defined by Eq. (|j). It amounts to the discontinuity of the momentum distribution 
at the Fermi surface and measures to content of correlations included in the model. As is well 
known the pairing modifies the chemical potential which is calculated self-consistently with the 
gap equation from the closure equation for the density of neutrons. In our approximation it is 
given by 

d 3 p r du 
(2tt) 3 J 2vr? 



P = 2 I 7^ I WZlG> + ) (14) 



d 3 p Z p u p 



( 2 -) 32vl v^+^ } (15) 

The prefactor 2 is due to spin degeneracy. Our numerical investigation is based on the solution 
of the two coupled gap equations, Eqs.(12) and (15), along with the self-energy approximated 
to the second order of the hole-line expansion. 

2. Results and conclusions. 

The Argonne AlA potential has been adopted as pairing interaction which is consistent with 
the self-energy data where the same force has been used. The gap equation has been solved 
in the form of Eq. (12), coupled with Eq. (15). This is a quite satisfactory approximation, 
especially in view of studying the self-energy effects. The results are reported in Fig. 3 for 
a set of different fci?-values. The solid line represents the solution of the gap equation in the 
standard BCS limit with the free single-particle spectrum. This is very close to the prediction 
obtained replacing the bare mass by the effective mass calculated in the BHF but still keeping 
Z = 1 (patterned line). This similarity stems from the fact that at the Fermi surface m*/m 
from BHF is close to one as shown in Fig. 2. The self-energy effects are estimated in two 
approximations. In the first one m* and the ^-factor are calculated from the approximation 
X = X 1 in a BHF code. In the considered density domain the Z-factor is around ~ 0.9. Despite 
its moderate reduction a strong suppression of the gap is obtained as shown in Fig. 3 (upper 
long-dashed line). It is due to the exponential dependence of the gap on all quantities. Still 
a further but more moderate reduction is obtained when the rearrangement term is included 
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in the second approximation E ~ X 1 + £ 2 (short-dashed line). The smaller ^-factor [Z ~ .83 
at kp = .8/m _1 ) is to a certain extent counterbalanced by an increase of the effective mass 
(m*/m m 1.2 at the same fcp). 

Self-energy corrections are mostly concentrated around the Fermi surface; therefore it is not 
appropriate to use Eq. (13) for a quantitative prediction of the gap, since it extends their effect 
beyond the Fermi surface. The self-energy effect turns out to be overestimated by Eq. (13) as 
we checked numerically. In Fig. 3 the results are reported for the two adopted approximations 
of the self-energy (lower long-dashed and short-dashed lines, respectively). 

In Ref. || the selfenergy effect has also been investigated within the generalized gap equation 
but the selfenergy has been considered only at the level of BHF approximation. However the 
approximations adopted in Ref. || for solving the gap equation understimate the correlation 
effects so that a moderate reduction is obtained. In the present approach the reduction is more 
pronounced with only S 1 . The inclusion of the rearrangement term S 2 brings about a further 
non negligible reduction. 

In conclusion we have shown that the superfluidity of a strongly correlated Fermi system 
requires to be described in the context of the generalized Gorkov approach. The fact that 
the quasiparticle strength can be significantly smaller than one cannot be counterbalanced by 
a corresponding enhancement of the effective mass. Moreover we have shown that reliable 
predictions from the generalize gap equation can only be obtained if the correlation effects 
are fully taken into account. We have treated here neutron matter but we expect that the 
selfenergy effects on the gap are very similar in symmetric nuclear matter. A next important 
step forward will be to include not only selfenergy effects but also, on an equal footing, vertex 
corrections. This shall be studied in a future work. 
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Figure captions 



• Fig.l One hole-line off-shell (solid line) and two- hole line off-shell (dashed 
line) at a fixed momentum k for three neutron densities p = kp/3n 2 . 

• Fig. 2 Effective mass (upper panels) and e-mass (lower panels) for three densities p = 
kp/3ii 2 in the BHF approximation (dashed line) and EBHF approximation (solid line). 

• Fig. 3 Energy gap in different approximations (A in the text). Solid line: free single- 
particle spectrum; patterned line: effective mass from BHF approximation and Z = 1; 
upper (lower) long-dashed line: solving the gap equation in the form of Eq. (12) (Eq. (13)) 
with BHF effective mass and Z; upper (lower) short-dashed line: solving the gap equation 
in the form of Eq. (12) (Eq. (13)) with EBHF effective mass and Z. 
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